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I. INTRODUCTION 

Ising models are paradigmatic in analytical and nu- 
merical studies of phase transitions [1, 2]. Their virtue 
is being simple enough to handle but nonetheless com- 
plex enough to capture the relevant physics. The same is 
true for other emblematic classical spin models, like the 
Potts model [3], or the six vertex model [4], which serve 
as toy models for certain physical systems. As a matter 
of fact, their applicability extends beyond physics, since 
spin models are used in the study of neural networks 
[5], biology [5] and, more generally, statistical mechani- 
cal tools are also applied in economics [6] . The reason lies 
in the fact that these models study classical degrees of 
freedom ( "spins" ) which interact with each other (possi- 
bly in many-body interactions), and this general scheme 
can serve as an abstract model for a wide class of systems. 

The central problem in the study of these models in 
equilibrium is the computation of their partition func- 
tion Z := e~^^^^\ where H{s) is the Hamiltonian 
(or: "energy function"), j3 is defined as {3 := 1/(/cbT) 
where /c^ is Boltzmann's constant and T is the tempera- 
ture, and s is the spin configuration. As is well known in 
statistical mechanics, this function captures all relevant 
physical properties, since any thermodynamical quantity 
(such as the magnetization or the mean energy) can be 
derived as a function of Z [7]. In other words, knowl- 
edge of >Z as a function of the parameters of the system 
amounts to complete knowledge of the thermal proper- 
ties of the system. Thus, it is natural to investigate the 
computational complexity (colloquially speaking: the ef- 
fort in terms of resources) required to compute or ap- 
proximate Z. The emblematic models mentioned above 
have received most of the attention in this direction. For 



example, for the Ising model, Barahona [4] showed that 
computing its partition function in three dimensions (3D) 
is computationally very hard — to be precise, it is #P- 
complete, which is, colloquially speaking, the counting 
version of NP [8]. This was an important contribution 
that settled the issue for those trying to tackle the prob- 
lem after Onsager's success in solving the Ising model in 
two dimensions (2D) [9]. 

Further, one can raise the question of how hard it is 
to compute the partition function of these models on a 
quantum computer. Quantum computers are known to 
offer a speedup over their classical counterparts in certain 
algorithms. Notably, the factoring problem is known to 
be feasible by a quantum computer [10], while it remains 
intractable in all known classical algorithms. Roughly 
speaking, 'feasible' means that the resources (in time and 
space) required to solve it scale polynomially with the 
size of the input of the problem, and intractable means 
that they may scale exponentially. The quantum com- 
putational complexity of classical spin models has been 
addressed, for example, in [11], where a quantum algo- 
rithm for the Ising partition function was presented (see 
also [12-16] for related work). Further, in [17] it was 
proven that computing the partition function of the Potts 
model is BQP-complete (see also [18, 19]). BQP stands 
for bounded-error quantum polynomial time, and, collo- 
quially speaking, is the class of decision problems that 
can be efficiently approximated by a quantum computer, 
and the hardest problems in this class are called BQP- 
complete. This situation contrasts with that of a different 
class of classical spin models, namely lattice gauge the- 
ories with gauge group Z2 [20], for which, to the best 
of our knowledge, no results are known concerning their 
quantum computational complexity. These are models 
with "Ising variables" (i.e. classical degrees of freedom 
with two states) but which nonetheless exhibit local sym- 
metries [21, 22]. 
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In this paper we tackle the question of the quantum 
computational complexity both of several paradigmatic 
classical spin models as well as of a Z2 lattice gauge the- 
ory. Our approach builds upon a mapping between par- 
tition functions and quantum circuits introduced in [23] . 
In that work, this mapping is exploited to show, among 
others, that estimating the partition function of the six 
vertex model and of Ising-type models is BQP-complete. 
Here we revisit this approach and extend the mapping to 
standard Ising models. Potts models and to Z2 lattice 
gauge theories. Based on that, we provide efficient quan- 
tum algorithms to estimate (with polynomial accuracy) 
the partition function of the six vertex model on a 2D 
square lattice, the Ising model with magnetic fields on 
a planar graph, the Potts model on a quasi 2D square 
lattice, and a 3D lattice gauge theory with gauge group 
Z2. Moreover, we show that computing the partition 
functions of these models is BQP-complete, that is, it 
is as hard as simulating arbitrary quantum computation. 
Therefore, in this work we put paradigmatic classical spin 
models and Z2 lattice gauge theories on an equal foot- 
ing as far as their quantum computational complexity is 
concerned. 

However, a word of caution is needed here: our re- 
sults are valid mostly for a complex parameter regime of 
the models. This means that we can prove the complex- 
ity results only if the coupling strengths of these models 
have certain imaginary values (this problem is also en- 
countered in [17, 23]). Note that although such complex 
parameters do not correspond to physical models, the 
partition function with complex arguments is commonly 
studied, e.g., in the context of evaluating the Tutte poly- 
nomial or finding (complex) zeros of Z to identify phase 
transition points [24]. We also want to stress that our 
results rely crucially on the fact that an additive approx- 
imation of the partition function is obtained (in contrast 
with a multiplicative approximation or an exact calcula- 
tion), which, roughly speaking, means that the desired 
quantity is approximated with polynomial accuracy. 

To summarize, the main results of this work are the 
following. 

Statement of results. Consider the partition func- 
tion Z of the following classical spin models 

(i) the six vertex model defined on a rectangular grid 

(a) the Ising model with magnetic fields defined on a 
planar graph, 

(Hi) the three-level Potts model on a quasi 2D square 
lattice with certain boundary conditions, and 

(iv) the Z2 lattice gauge theory on a 3D square lattice, 

defined on a certain complex parameter regime, that is, 
the value of coupling strength J is, for example, e^^ = i. 
Furthermore, consider that there are certain values of the 
coupling strengths which appear together (i.e. certain sets 
of neighboring spins whose interaction takes a specific 
value). Then we provide efficient quantum algorithms to 



approximate the partition function Z of these models with 
polynomial accuracy. Moreover, we show that estimating 
these partition functions is BQP-complete, that is, it is 
as hard as simulating arbitrary quantum computation. 

Finally, as an extension of our results, we show that, if 
the Ising model considered above is defined on a lattice 
with periodic boundary conditions, then estimating its 
partition function with polynomial accuracy is DQCl- 
hard. This is the representative class of a scheme for 
quantum computation called the 'one clean qubit model', 
where all qubits but the first one are initialized in a to- 
tally mixed state [25]. 

This paper is structured as follows. In Sec. II we give a 
background on classical spin models, as well as on some 
notions of complexity of classical spin models. Then we 
review how to define BQP-complete problems as an esti- 
mation of a unitary matrix element (Sec. III). In Sec. IV 
we show how to relate partition functions of classical spin 
models with spin circuits, based on [23] and extended 
here. The main results of this work are presented in 
Sec. V, where we show the BQP-completeness of the six 
vertex model, the Ising model, the Potts model, and of 
the 3D lattice gauge theory with gauge group Z2 with 
certain conditions. In Sec. VI we present an extension of 
the results related to the one-clean-qubit model. Finally, 
we present our conclusions in Sec. VII. 

II. CLASSICAL SPIN MODELS 

In this section we will present some general consider- 
ations about classical spin models as well as some facts 
concerning their computational complexity. 

A. General background 

Classical spin models have proven successful to model 
magnetism, where they capture interesting physics such 
as critical phenomena despite their simplicity. More gen- 
erally, classical spin models can serve as toy models for 
complex systems. For example, Ising models have been 
used to model neural networks (see Hopfield networks), 
or spin glasses (see, e.g., [5]). 

Let us now define what is understood by a classical 
spin model. Common to all of them are the following 
ingredients: 

(i) A degree of freedom represented by a classical spin 
which may take on a set of values: s G {0, 1, g — 
1}. This is a g'-level system. 

(ii) A lattice, or more generally an arbitrary graph G, 
to which the classical spins are associated. They 
can sit at the vertices, edges or faces of the graph; 
the idea is that the graph encodes the interaction 
pattern of the model. The complexity of the model 
is affected by the lattice or graph on which it is 
defined, as we will see below. 
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(iii) An energy function H{s) depending on a given spin 
configuration and the coupling strengths represent- 
ing the types of interactions of the model: nearest- 
neighbor, many-body interactions, magnetic fields, 
etc. 

(iv) A partition function Z which is obtained by sum- 
ming over the Boltzmann weights of all spin config- 
urations, Z = ^-^H{s) ^ 

From this common structure, several different families 
of models can be distinguished. One of the most relevant 
criteria is whether the model exhibits global or local sym- 
metries (if any). We shall refer to the former as standard 
statistical models and to the latter as lattice gauge the- 
ories. This distinction of symmetry has profound conse- 
quences in the physics of the models. It is also related to 
their range of applicability: standard statistical models 
appear naturally in descriptions of classes of condensed 
matter systems and lattice gauge theories originated in 
the study of the fundamental interactions in nature and 
elementary particles. 

Dimensionality is another key distinction which deter- 
mines the complexity of the models. However, note that 
in order to have a well-defined notion of dimension for a 
graph G, this must be embedded in a smooth manifold 
of dimension D. 

Standard statistical models can be divided into two 
big families depending on where the interactions on the 
graph G take place: vertex models and edge models. Ver- 
tex models were introduced to describe ice-type models, 
crystals with hydrogen bonding or ferroelectrics [4] . Edge 
models were introduced to explain phase transitions in 
materials with elementary magnetic moments [26, 27]. 

A vertex model consists of classical spins, namely q- 
level particles Se G {0,1,...,^' — 1}, which are placed on 
the edges of the lattice, and (typically many-body) in- 
teractions take place on the vertices. In the case of a 
(tilted) 2D lattice, one deals with four-body interactions 
between neighboring particles, and each of the spins par- 
ticipates in only two interactions. The Hamiltonian of 
such a system is given by 

H =^ h''{si,Sj,Sk,si) , (1) 

where h^{si^ Sj^ Sk-, si) is a (local) four-body interaction 
term between the spins Si^Sj^s^^si sitting at the edges 
incident to vertex a. We denote the Boltzmann weight 
associated with this local energy as 

w^{si,sj,sk,si) := e-'''^"(^-^-«'='^') . (2) 

The partition function is obtained by multiplying all lo- 
cal Boltzmann weights, and summing these over all spin 
configurations s, 

^™ = ^ JJ^^;''(s^,s^•,5/e,s^) . (3) 

s a 



In contrast to vertex models, in edge models the clas- 
sical spins sit at the vertices of the graph, also taking q 
possible states, Si G {0, . . . , (7 — 1}, and interactions take 
place along the edges. Consider, thus, a g-state edge 
model on an n X m square lattice with an edge-dependent 
energy function h^{si,Sj). Let 

w%s,,Sj):=e-^''^^'-'^^ (4) 

denote the corresponding Boltzmann weight. Then the 
partition function is given hy Z = Ile^ij '^^i^i^ ^j)- 

Concerning Lattice Gauge Theories (LGTs), the most 
relevant criterion to classify them is whether the inter- 
nal gauge group Q is Abelian (discrete like Q = ^g, or 
continuous Q = U{1))^ or non- Abelian (discrete such 
as a permutation group Q = 83^ or continuous such as 
Q = SU{N)). We refer the reader to [20] for an in- 
troduction to these models. In this work we will fo- 
cus on Z2 LGTs, and we will be in the following fea- 
tures: they are models whose classical spins can take 
two values Se G {0, 1}, they sit at the edges of a d- 
dimensional square lattice, and they interact along the 
faces of this lattice. More precisely, the interaction of 
spins Sj^ si at the boundary of face f^df^ has the 
form 

hf{Si,Sj,Sk,Si) = -JfS{Si^Sj +S/e+S^), (5) 

where the sums are performed modulo 2 throughout this 
section, and 5{0) = 1 and it is otherwise. The Hamilto- 
nian is then obtained as a sum over interactions on every 
face: 

H{s) = -Y,hf{{se:eedf}), (6) 
/ 

where df denotes the boundary of face /. 

These models exhibit Z2 gauge symmetry; more pre- 
cisely, its Hamiltonian is invariant under Z2 operations 
around any vertex, gy = neGinc^;^e, where e e inc v 
denotes all edges incident to vertex and Xg is a fiip 
operator, Xg : s ^ 5 + 1. One can use this symmetry to 
eliminate some degrees of freedom, a process usually re- 
ferred to as 'gauge fixing'. A specific choice of this fixing 
is the 'temporal gauge', where all degrees of freedom in 
one particular direction (the one associated with time) 
are fixed. A restriction about gauge fixing that concerns 
us is the fact the edges whose variable has been fixed by 
the gauge cannot form a closed loop [28]. This fact will 
be important in our proof of the BQP-completeness of 
this model in Sec. VD. 



B. Classical computational complexity of spin 
systems 

In this section we will be interested in the classical 
computational complexity of the classical spin models 
presented in Sec. II A. Generally speaking, understanding 
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the properties of, say, the Ising model on some graph is a 
difficult task. This is reflected by the fact that the Ising 
and other models are associated with hard problems in 
computational complexity theory [4]. For concreteness 
we will focus in the following on the Ising model, but the 
considerations in this section have general relevance. 

Most prominently, the Ising model is known to be 
associated with computational problems which are NP- 
complete; here NP stands for 'non-deterministic polyno- 
mial time'. The complexity class NP consists of all deci- 
sion problems / (i.e. YES/NO questions) which have the 
property that, for every input x for which it is claimed 
that f{x) = YES, there exists a 'short proof that this 
is indeed the case, i.e. a proof which may be efficiently 
verified as a function of the size of the input. More pre- 
cisely, for each problem / G NP it is required that there 
exists an efficiently computable function V (the verifier 
of the proof) such that: 

For every input x, one has f{x) = 'YES' if and only 
if there exists a poly-size bit string ^ (the witness) 
satisfying V(x,0 = 'YES'. 

Colloquially speaking, NP problems have the property 
that, whenever a solution to the problem is proposed 
(e.g. by an untrusted third party), it is possible to ef- 
ficiently verify whether this proposed solution is indeed 
correct. Note that in the definition of NP no mention- 
ing is made of the difficulty of finding a solution; even 
though a problem has an efficient verifier, it is a priori 
not excluded that the time required to find a solution 
scales exponentially with the input size. 

An archetypical N P problem related to the Ising model 
is the problem of deciding whether the ground state en- 
ergy of Hg{s) on a graph G (which constitutes the input 
of the problem) is below a certain value K. This prob- 
lem is indeed in NP: if the ground state energy of H{s) 
is smaller than then the ground state provides a wit- 
ness which allows to efficiently verify this fact; the verifier 
function is nothing but the energy function Hg{s). 

Not only is the problem of determining the Ising 
ground state in NP, it is among the hardest problems in 
this complexity class. This is refiected by the fact that 
this problem is known to be NP-complete. This means 
that every problem in the class NP can be reduced, with 
only polynomial computational effort in the input size 
of the problem, to an instance of the Ising ground state 
problem. This implies, in particular, that the existence 
of an efficient algorithm for the Ising ground state prob- 
lem would yield an efficient algorithm for all problems in 
NP. The N P-completeness of the Ising model thus points 
to an intrinsic computational difficulty of this simple sys- 
tem. 

There are several variants of the Ising ground state 
problem which are known to be NP-complete. We men- 
tion two of them. 

Theorem 1 [4] The following problems are NP- 
complete: 



(i) Given a graph G and an integer K, deter- 
mine whether the ground state energy of H{s) = 
Xle=o6 ^a^b ^-^ Smaller than K. 

(a) Given a planar [40] graph G and an integer K, de- 
termine whether the ground state energy of H{s) = 
^e=ab ^aSb + Z)a ^cl Smaller than K. 

We remark that, in the second of these problems, the 
presence of the external fields {ha = —1) is crucial to 
obtain N P-completeness. Indeed, it is known that the 
ground state energy, as well as the partition function, 
of the Ising model on an arbitrary planar graph with- 
out external fields can be efficiently computed. We also 
note that the quantum computational complexity of the 
Ising model, as will be discussed in Sec. VB, will involve 
Ising models on planar graphs in the presence of magnetic 
fields. 

The N P-completeness of the above ground state prob- 
lems has strong implications for the evaluation of the 
corresponding partition functions. First, once the parti- 
tion function of a model can be evaluated efficiently, also 
the ground state energy of the model can be efficiently 
determined: the evaluation of the partition function is 
'at least as hard' as the evaluation of the ground state 
energy. Consequently, for the NP-complete Ising mod- 
els, the evaluation of their partition function is NP-hard, 
i.e. at least as hard as any problem in NP. Note, how- 
ever, that the evaluation of the partition function does 
not belong to the class NP, as it is not a decision problem 
but rather a counting problem. The relevant complexity 
class in this case is #P (' sharp- P'). Given an efficiently 
computable decision problem g (i.e. ^ G P), the problem 
of determining how many inputs yields the answer 'YES', 
represented by the number ^g = \{x : g{x) = 'YES'}|, 
defines the complexity class #P. 

Since the ground state problems of the Ising models 
in Theorem 1 are NP-complete, it can be shown that 
computing the corresponding partition functions are #P- 
complete problems: every problem in #P can be reduced, 
with polynomial computational effort, to the evaluation 
of the partition function of such an Ising model on some 
graph. This is formulated in the following result: 

Theorem 2 [4] The following problems are #P- 
complete: 

(i) Given a graph G and A = , determine the par- 
tition function Z{X) of the Ising model on G with 
energy H{s) = Y.e=ab^aSb' 

(a) Given a planar graph G and A = ^ determine 
the partition function Z{X) of the Ising model on 

G with H{s) = Y^e^ab ^aSb + Ea "^a- 

III. UNITARY MATRIX ELEMENTS AND 
BQP-COMPLETENESS 

The goal of this paper is to relate the evaluation of par- 
tition functions to problems which are complete for quan- 
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turn complexity classes. The main complexity class which 
will be considered is 'bounded-error quantum polynomial 
time' (BQP), representing the class of decision problems 
which can be solved efficiently on a quantum computer. 
The route we will take to prove BQP-competeness of cer- 
tain partition function problems will be to start from a 
standard complete problem for BQP and then to relate 
these problems to the approximation of partition func- 
tions. The standard BQP-complete problem in question 
involves estimating matrix elements of unitary quantum 
circuits, which we briefly discuss here. 

Consider a quantum circuit U acting on n qubits, 
which is composed of poly(n) gates acting each on, say, 
at most two qubits. Let {|0), |1)} denote the single-qubit 
computational basis. Then there exists a well known 
technique to estimate the matrix element (0|*^^[/|0)*^"^ in 
poly-time on a quantum computer, using the 'Hadamard 
test' (see Fig. 1 and its caption for a brief a explanation — 
we refer to the literature, e.g. [18, 19, 29, 30], for further 
explanations). More precisely, for any approximation 
scale e which scales at most inverse polynomially with 
n, the Hadamard test returns a complex number c which 
satisfies 



(7) 



with a success probability that is exponentially (in n) 
close to 1. 
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FIG. 1: The Hadamard test. This is performed with a 
quantum circuit where the first qubit is transformed with a 
Hadamard gate, then the unitary U is applied to rest of the 
qubits conditional on the state of the first qubit, and then a 
Hadamard gate is applied to the first qubit again. Finally, 
this qubit is measured in the az basis. From the probabil- 
ity to obtain and 1, po and pi respectively, one can esti- 
mate the real part of c (denoted Re(c)), as they are related 
by Po = [1 + Re(c)]/2 and pi = [1 - Re(c)]/2. To estimate 
the complex part, the circuit is modified to include a phase 
gate P = diag(l,i) after the first Hadamard gate and be- 
fore the controlled-^/ gate. The probabilities to measure 
and 1 are then related to the imaginary part of c (denoted 
Im(c)) as Po = [1 — Im(c)]/2 and pi = [1 + Im(c)]/2 (see also 
[18, 19, 29, 30]). 



Moreover, estimating the above matrix element prob- 
lem is BQP-hard, i.e. every decision problem that can 
be solved efficiently with a quantum computer can be 
reduced, with polynomial (classical) computational ef- 
fort, to the estimation of such a matrix element with the 
aforementioned accuracy e. Without loss of generality we 



can restrict all operations to act on nearest-neighboring 
qubits, since the SWAP operation can be used to move 
distant qubits to nearest-neighbor positions (with linear 
overhead in the number of qubits). Thus, one has the 
following. 

Theorem 3 The following problem is BQP-hard: 

Consider an n-qubit quantum circuit U consisting 
of a polynomial number of two-qubit gates acting on 
nearest-neighbor qubits. Then provide a number c 
such that 



(0|^^/7|0) 



(8)n I 



< 



poly(n) 



(8) 



holds with a probability which is exponentially (in 
n) close to 1. 

We provide a few remarks: 

(i) One may adapt the formulation of Theorem 3 in a 
straightforward way to arrive at a BQP-complete 
decision problem (i.e. a problem that is both BQP- 
hard and in BQP). To do so, one considers circuits U 
for which it is promised that |(0|^^/7|0)^^| is either 
< 1/3 or > 2/3, and the goal is to decide which of 
these two cases holds. 

(ii) BQP-hardness of the matrix element problem in 
Theorem 3 is maintained if one considers quantities 
of the form (?/;|[/|?/^'), where and are fixed 
complete product states, instead of {{)\^'^U\G)®^ . 
This is because instances of the former matrix ele- 
ment problem can easily be reduced to instances of 
the latter. This property will be used in Sees. VA, 
VB and VC. 

(iii) BQP-hardness is maintained when the circuit U is 
restricted to be composed of gates from a strongly 
universal elementary gate set. An elementary gate 
S consisting of gates acting on at most two qubits is 
said to be strongly universal if every two-qubit uni- 
tary operation can be approximated with arbitrary 
accuracy by products of elements from S. Further, 
BQP-hardness is also maintained if, instead of con- 
sidering strongly universal unitary gate sets, gate 
sets are considered that are encoded universal for 
quantum computation, for a suitable notion of en- 
coded universality. This will be important for our 
proofs in Sees. VA, VC and VD. 



IV. MAPPINGS BETWEEN CLASSICAL 
LATTICE MODELS AND QUANTUM CIRCUITS 

In this section we review the mappings between classi- 
cal lattice models and quantum circuits (introduced and 
used in [23]), and we will extend them to standard Ising 
models, to Potts models and to Z2 LGTs. The map- 
pings will allow us to interpret the partition function 
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of the classical model as a quantum expectation value. 
More precisely, consider a circuit C consisting of (uni- 
tary) quantum gates. Then we will show that partition 
functions Z of vertex models, edge models and Z2 LGTs 
are related to matrix elements of certain quantum cir- 
cuits C as 



i^{L\C\R). 



(9) 



where {L\ {\R)) are product states determined by the left 
(right) boundary conditions of the classical spin model 
(see below), and /t: is a constant depending on the size of 
the lattice, i.e. the number of classical spins. 



A. Vertex models 

We start by considering vertex models (see Sec. II A). 
Here we essentially follow the argument of [23]. For il- 
lustration purposes, we will concentrate on a tilted 2D 
square lattice. However, our mappings are not restricted 
to such lattices and are easily generalized to other (reg- 
ular) lattices. 

Our construction begins with the following observa- 
tion. The local Boltzmann weights (2) associated to each 
interaction can be seen as rank-four tensors of dimension 
g', or, equivalently, as q'^ x q'^ matrices by grouping the 
indices into left indices (i^j) and right indices (/c,/). We 
take the latter approach to associate each of these matri- 
ces with a quantum gate acting on two g'-level quantum 
states (see Fig. 2), 



Si,Sj,Sk,Si 



W^'iSi, Sj.Sk, Si)\Si, Sj){Sk, Si I . (10) 



Note that the right indices (k^l) correspond to the input 
of the gate, while the left indices (i^j) represent the out- 
put. The state is thus processed from right to left, where 
gates corresponding to outermost right vertices are per- 
formed in the first place [41]. The corresponding quan- 
tum circuit C is given by m layers of nearest-neighbor 
two-qubit gates (see Fig. 2), 



(11) 



This can be seen as the contraction of a tensor network, 
i.e. as a summation over joined indices. For translational 
invariant models, each of the layers corresponds to the 
transfer matrix [4] of the classical model. 

Thus, generally speaking, we have mapped a product 
of interactions of the classical spin model (which is es- 
sentially a partition function) to a contraction of quan- 
tum gates (which is essentially a quantum circuit). Now 
we only need to show how to map the boundary spins. 
We consider an n x m lattice with fixed boundary con- 
ditions, namely a lattice whose spins at the left and 
right boundary are fixed in some arbitrary configuration 



L = (sf , . . . , s^) and R = (sf , . . . , 5^), respectively. Us- 
ing the spin states as computational basis states, we map 
these to two n-particle quantum states: 

{L\ = {s^\...{sfi\, 
\R) = |sf)...|s«). 



(12) 



where we omit the tensor product symbol throughout 
this paper. Note that, since the qubits are processed 
from right to left, the state \R) serves as input for the 
circuit, while \L) constitutes the readout basis state. 

It is now straightforward to see that the overlap of the 
resulting state C\R) with the product state {L\ is exactly 
the partition function of the classical vertex model. 



{L\C\R). 



(13) 



This concludes the mapping between matrix elements of 
quantum circuits and partition functions of classical spin 
models. 
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FIG. 2: (Left) In a vertex model, particles (black dots) sit 
at the edges and interactions (pale red dots) take place in the 
vertices. (Right) This model is mapped to a quantum circuit, 
where each particle becomes a qubit, and each interaction a 
two-qubit gate. 

Now we make some remarks concerning this mapping. 

(i) Complex couplings. Note that the quantum gates 
are specified by the parameters of the classical 
model, namely by the Boltzmann weights of the 
local interactions. It follows that the gates 
(Eq. (10)) are unitary only in certain parameter 
regimes of the interaction of the classical spin 
model. This will lead to the requirement of com- 
plex parameters in our proof of Sec. V A. 

(ii) Open and periodic boundary conditions. This map- 
ping can be easily extended to open boundary con- 
ditions, i.e. to systems where the left and right spins 
are 'free' and thus fully summed out in the partition 
function. This is achieved by replacing the left and 
right states \L) and \R) by the state |+)*^^, where 
1+) = q~^^'^ '}Zi=o\^) ^ superposition over all q 
single spin states. This gives rise to the identity 

Z^^C =q"{+\^^C\+)'»". (14) 

Furthermore, periodic boundary conditions can also 
be taken into account by summing over the diagonal 
matrix elements, which results in 

ZZ^=Tv{C). (15) 
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(iii) Other geometries. One can also consider other ge- 
ometries such as vertex models on a 2D tilted trian- 
gular lattice [31], where six-body interactions take 
place at the vertices and the Boltzmann weights can 
be arranged into x matrices, corresponding to 
a quantum circuit with three-body quantum gates. 
One such model is the 32 vertex model [31]. Also 
3D models, such as models on a tilted 3D square 
lattice, are of this type and can thus be mapped 
to quantum circuits. In this case, one deals with 
three-body gates acting on a 2D array of quantum 
particles. 



B. Edge models 

In the following we will present a mapping for edge 
models, also following [23]. Unlike vertex models, in 
edges models interactions take place at the edges, as ex- 
plained in Sec. II A. For the mapping we will distinguish 
between interactions at horizontal or vertical edges. More 
precisely, we associate di q x q matrix to each horizontal 
edge e: 



(16) 



and 8i q^ X q^ diagonal matrix to each vertical edge e. 



'= '^W%Si,Sj)\Si,Sj){Si,Sj 



(17) 



The matrices Wl^ and will be regarded as (possi- 
bly non-unitary) quantum gates acting on a single, re- 
spectively a pair, of g-level quantum systems. We now 
consider a ID quantum system composed of n g-level 
systems and the quantum circuit C acting on this system 
as depicted in Fig. 3. The circuit C consists of alternat- 
ing layers of operations associated with the horizontal 
and vertical edges of the 2D lattice. Each round of C 
associated with a layer of horizontal edges consists of a 
product of one-local operators VFg^, whereas every round 
associated with a layer of vertical edges is a product of 
(commuting) two-local operations . We define com- 
putational basis states \L) and \R) associated to the left 
and right boundary conditions, respectively, analogously 
to Eq. (12). With these definitions, one has the following 
correspondence: 



rL,R 



{L\C\R) . 



(18) 



This equation is readily verified by employing the defini- 
tions of gates and . 

We emphasize that the comments concerning the map- 
ping for vertex models apply to this mapping with 
straightforward modifications: unitary gates lead to com- 
plex coupling strengths (see Sec. IV A remark (i)), and 
the mapping can be extended to open and periodic 
boundary conditions (see Sec. IV A remark (ii)). Now 
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FIG. 3: (Left) In an edge model, particles (black dots) sit at 
the vertices and interactions (pale red and blue ellipses) take 
place along the edges. (Right) This model is mapped to a 
quantum circuit, where particles are mapped to qubits, inter- 
actions along the time direction become single-qubit gates, 
and interactions perpendicular to time become diagonal two- 
qubit gates. 



we give some further modifications of the mapping for 
edge models that will be useful for the following sections. 

(i) Consider that at site i in the lattice a local magnetic 
field is present. This is represented by an additional 
term hi{si) in the energy function and correspond- 
ing Boltzmann weight 



(19) 



with Si = 0,...,(7 — 1. We associate to it the follow- 
ing diagonal q x q matrix 



Wi := ^Wi{si)\si){si\ . 



(20) 



Now a mapping to a quantum circuit C can be estab- 
lished in a similar fashion as above, with the distinc- 
tion that each layer associated with a slice of verti- 
cal edges now consists of a product of the associated 
two-qubit gates and the associated single-qubit 
gates Wi. Note that, as all such gates are diagonal 
operations, there is no problem regarding operator 
ordering. With this choice of C, it can readily be 
verified that the associated partition function can 
be written as (18). 

(ii) These mappings may be easily generalized to graphs 
other than the 2D lattice, also similarly as for vertex 
models (see Sec. IV A remark (iii)). In particular, 
below we will consider the following class of sub- 
graphs of the 2D square lattice: a graph G is said to 
be a planar circuit graph if it can be obtained from 
an n X m rectangular grid (for some n and m) by 
deleting a subset of vertical edges and contracting 
a subset of horizontal edges. We call n the vertical 
dimension of G; note that this quantity is uniquely 
defined for every planar circuit graph. Similar as for 
the 2D square lattice, one can associate a quantum 
circuit C with every planar circuit graph; more pre- 
cisely, such circuit acts n g'-level systems, and one 
associates each horizontal and vertical edge e with 
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the gates and ^ respectively. Furthermore, 
local magnetic fields acting on the particles can also 
be easily incorporated by associating a gate Wi with 
each vertex z; see Fig. 4 for an example. This map- 
ping will be relevant in Sec. VB. 



Wi 
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FIG. 4: (Left) Edge model with magnetic fields defined on a 
planar circuit graph. The latter is obtained from a rectangu- 
lar grid by deleting some vertical edges and contracting some 
horizontal ones. (Right) This model is mapped to a quantum 
circuit: particles are mapped to qubits, and horizontal and 
vertical edge interactions are mapped to single-qubit (non- 
diagonal) and two-qubit diagonal gates, respectively, and lo- 
cal interactions (e.g. magnetic fields) are mapped to single- 
qubit diagonal gates. 

In this paper we will be interested in two particular 
edge models: the Ising model and the Potts model. We 
now specialize the above discussion to the case of the 2D 
Ising model. We consider the Ising with magnetic fields 
defined on (sublattices of) a 2D square lattice. The inter- 
action between spins Si and Sj located at the endpoints 
of edge e = (i, j) is given by 

he{Si, Sj) = -JeS{Si + Sj) , (21) 

and the contribution of the magnetic field at site i is 

hi{si) = -hiS{si). (22) 

Here the spin states Si may take values and 1, the sums 
are performed modulo 2 and as before ^(0) = 1 and it 
is otherwise. Further, Jg and hi are constants which 
represent the strengths of the pairwise interaction and 
magnetic field, respectively. With the definitions (16), 
(17) and (20), we have 



1 

1 



1 



WI = diag(e^^% 1,1,6^-^-), 



(23) 



Now we concentrate on the mapping for the Potts 
model. Given a graph G with vertex set V and edge 
set the Potts model [3] consists of g'-level particles 
sitting at the vertices of G and interacting along the 
edges of G. Let v denote two such g-level particles 
u^v G {0,1,... ,(7 — 1}, which interact along edge e. Then 
the Potts-type interaction is of the form 



where 6{0) = 1 and it is otherwise. We shall later 
consider a more general form of this interaction: 

h^'iu, v) = -Ju^yS{u -v) - Ju^v{l - S{u - v)) . (25) 

This amounts to a shift of the interaction energy which 
does not change the physics. Let s denote the spin state 
of all particles: s = (ii, . . .). Then the Hamiltonian of 
the Potts model is a sum of these two-local terms over 
all edges 



H{s) 



eeE 



h^{u, v) . 



(26) 



h^{u^ v) = —JeS{u — v) . 



C. Lattice gauge theories 

Now we focus on another family of models, namely Z2 
LGTs (see Sec. II A), and we will introduce mappings for 
their partition functions. 

We shall restrict the following discussion to a 3D Z2 
LGT with the temporal gauge. As explained in Sec. II A, 
fixing this gauge is achieved by fixing all spins lying on 
edges with a specific direction of the lattice (the 'time' 
direction). However, the mapping can be generalized in 
a straightforward manner to 3D Z2 LGTs with another 
gauge fixing of the spins — we will return to this com- 
ment below and in Sec. VD. Because of the temporal 
gauge fixing, interactions in the time direction are two- 
body interactions, whereas those in the spatial direction 
(i.e. in faces without an edge in the time direction) re- 
main four-body interactions as originally. To construct 
the mapping, we proceed similarly as above. We asso- 
ciate the Boltzmann weight of the two-body interaction 
at the temporal face 

w}{s,,s,):=e-^''f^^-^^\ (27) 

with a single-qubit (non-diagonal) gate Wj, 

Further, the Boltzmann weight of the four-body interac- 
tion at the spatial face 

w}{s„s„Sk,si) := e-'''^/(^-^-^-«') (29) 

is mapped to a four-qubit diagonal gate 

•= ^ wj{Si,Sj,Sk,Si) X 

Si,Sj ,Sk,Si 

\si,Sj,Sk,si){si,Sj,Sk,si\ . (30) 

Thus, this maps the 3D Z2 LGT to a quantum circuit 
where a 2D array of qubits is processed in the time direc- 
tion with single-qubit gates, and in the spatial direction 
with (diagonal) four-qubit gates (see Fig. 5). As before, 



(24) it follows that the partition function Z = J2s 11/ ' 
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is mapped to a quantum a circuit C = H/ ^/^/- 

L, R denote the left and right boundaries, respectively, 

as in Eq. (12). Then our mapping reads 



Z^^^ = {L\C\R) . 



(31) 



Similarly as for vertex and edge models, we note that 
unitary gates will be translated to complex coupling 
strengths (as in Sec. IV A remark (i)), and one can ob- 
tain similar mappings for LGTs with open and periodic 
boundary conditions (see Sec. IV A remark (ii)). Next we 
make some further comments about this construction. 

(i) The mapping can easily be extended to Z2 LGTs 
with other gauge fixings. As a matter of fact, we 
will use one of these mappings in Sec. VD. Since 
one deals with square lattices, one can always as- 
sociate one dimension of the lattice to time. Then, 
generally, one fixes some spins in the time direction 
and some others in the space direction, the only 
restriction being the avoidance of closed loops (see 
the comment in Sec. II A and [28]). This introduces 
new type of interactions. For example, there can be 
a face in the space direction where three spins have 
been fixed by the gauge, and thus only one variable 
is left: 



(32) 



Such interactions are mapped to single-qubit diag- 
onal gate W^: 



Wf = ^«;«(si)|si)(si|. 



(33) 



Examples of these gates will be given in Sec. VD. 
Note, however, that the mapping does not need to 
be defined on all faces: for the BQP-completeness 
proof it suffices to specify the coupling strength on 
every face. In particular, we will set J = on faces 
in the time direction where the temporal gauge is 
not fixed, and the interactions are not mapped to a 
gate in this case. 

(ii) Notice that, due to the form of the interaction of 
this system (given by (5)), each gate is specified by 
only one parameter: e^^. We will make repeated 
use of this fact in Sec. VD. 

(iii) This mapping can be extended to Zq LGTs (see [20] 
for an introduction of these models). In this case, 
each interaction can take q values, and they would 
translate to one- and foui-qudit gates, where each 
qudit is a g-level quantum system. 

(iv) More generally, the mapping can also be extended 
to Z2 LGTs defined on square lattices in d dimen- 
sions. The constructions is a generalization of the 
3D case, where one has to select a translational 
invariant time direction out of the d dimensions. 



and fix the temporal gauge. Since interactions also 
take place on the faces, time and space interactions 
are also mapped to one- and four-qubit gates, re- 
spectively. This results in a quantum circuit that 
processes a (d — 1) dimensional array of spins with 
single-qubit gates in the time direction, and with 
four-qubit gates in the space dimensions. Similar 
considerations apply to Zq LGTs, with one- and 
four-qudit gates. 



V. BQP-COMPLETENESS RESULTS 

This section contains the central results of this paper: 
here we will provide efficient quantum algorithms to es- 
timate the partition function of the six vertex model, the 
Ising model, the Potts model and the 3D Z2 LGT in a 
certain (complex) parameter regime, with polynomial ac- 
curacy. Moreover, we will show that approximating these 
partition functions is BQP-complete. 



A. Six vertex model 

Our goal is to prove that approximating the partition 
function Z of some vertex models in a certain (complex) 
parameter regimes is BQP-complete (see also [23]). To 
show that, we will make use of the mappings between ver- 
tex models and quantum circuits described in Sec. IV A. 

We consider the q = 2 six vertex model (or: 'ice-type 
model') and the eight vertex model [4] on a (tilted) 2D 
square lattice. In the six vertex model, only 6 of the 
16 possible spin configurations give rise to a non-zero 
Boltzmann weight. More precisely, is a 4 x 4 matrix 
of the form 



W = 



( '^00,00 




V 










^01,01 ^01,10 
^10,01 




^10,10 






(34) 



where we use 



si as a shorthand notation for 



w{si^ Sj^ Sk^ si) (see Eq. (10)). The eight vertex model 
[4] is obtained by additionally allowing the entries 
'^00,117 '^11,00 to be non-zero. We consider a parame- 
ter regime of the classical model where all matrices 
are unitary. This gives rise to a unitary circuit C formed 
of two-qubit quantum gates. Notice that this generally 
corresponds to (non-physical) complex parameters for ei- 
ther coupling strengths J or the inverse temperature /3 
as we pointed out in Sec. IV A remark (i). Finally, we 
assume that we have staggered left an right boundary 
conditions of the form L = R = (0101 . . . ). Our result is 
the following. 

Result 1 (BQP— completeness of the six vertex 
model) Consider the six vertex model defined on a 
n X polyin) rectangular grid with fixed boundary condi- 
tions. Further, consider that this is defined at inverse 
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FIG. 5: (Left) In a 3D LGT, particles (black and gray dots) sit at the edges and interactions (pale red and blue ellipses) take 
place on the faces. Gray dots indicate particles whose state has been fixed by the gauge. (Right) This model is mapped to a 
quantum circuit, where interactions along the time direction become single-qubit gates (blue squares), and those perpendicular 
to it become four-qubit gates (pale red squares). Thus, a 2D array of qubits is processed by a circuit consisting of a sequence 
of single-qubit gates and diagonal four-qubit gates. 



temperature /3 and with couplings strengths 



'^00,00 
'^^01,01 
'^^01, 10 



^11,11 

^10,10 
^10,01 



J2t 



cos(2t) 
i sin(2t) , 



(35) 



where t is a continuous parameter, and 



'^00,00 
'^^01,01 



^11,11 

^10,10 



1 

^01,10 



-^10,01 



1 

V2' 



(36) 



Let Z denote the partition function of this model. Then 
we provide efficient quantum algorithms to estimate Z 
with polynomial accuracy. We also show that the problem 
of approximating Z is BQP-complete. 

Before starting the proof, we remark boldface symbols 
will denote encoded states and operators throughout this 
paper. To prove Result 1, we will show that any quan- 
tum computation can be reduced to the evaluation of 
the partition function of a six vertex model on a tilted 
2D square lattice with staggered boundary conditions. 
We prove this statement in the following steps: 

1. We show that quantum gates of the form (34) are 
computational universal for encoded quantum com- 
putation. To do so, we use the four-qubit encoding 
for |0) given by Refs. [32, 33], which is of the form: 



io) = i(ioi)-iio)r. 



(37) 



Note that |1) can be prepared by means of the en- 
coded universal circuit that we will show next. 

Now we consider the exchange (or Heisenberg) in- 
teraction, 

^ex = CTx ^ ^ ay ^ O-y ^ ^ , (38) 

with corresponding two-qubit gates 



The Heisenberg interaction is (encoded) universal 
for quantum computation [32, 33]. In other words, 
by using gates of the form (39), one can prepare 
any quantum state = C|0) in an encoded form. 
Gates of the form (39) can be generated with six 
vertex-type gates, i.e. of the form (34), by setting 
the non-zero entries specified in (35). 

2. Now we show that the encoded initial state |0)*^^ 
can be prepared from the state corresponding to 
staggered boundary conditions |0101...). To this 
aim, we consider an operation V of the form (34) 
with the non-zero entries of Eq. (36). It is straight- 
forward to check that 1/|01) = (|01)-|10))/V2, and 
hence 



0) = v^^ioioi). 



(40) 



JtHe 



(39) 



3. Finally, we observe that matrix elements of the 
form (0|^^C|0)^^ can be efficiently approxi- 
mated by a quantum computer with polynomial ac- 
curacy, as long as C can be implemented efficiently. 
Since we are dealing with a poly-size quantum cir- 
cuit consisting of two-qubit gates (as indicated in 
Eq. (34)), C can be implemented efficiently. The 
estimation of the matrix element is achieved by the 
Hadamard test, as pointed out in Sec. III. These 
overlaps are related to partition functions of the 
six vertex model in a certain (complex) parameter 
regime via (13), since all gates [/, V involved are of 
the form (34), and thus correspond to Boltzmann 
weights of the six vertex model. 

This concludes the proof of Result 1. 

We make a few remarks concerning this construction. 

(i) Note that the complex parameter regime of 
Eqs. (35) and (36) is due to the fact that these 
entries correspond to unitary gates, as noted in 
Sec. IV A remark (i). 

It is worth mentioning that in [34] it is shown that 
universal quantum computation can be achieved 
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with real gates alone (i.e. gates with real entries). 
Thus, at first sight, it may seem that applying this 
method to our circuits would allow us to prove re- 
sults for a real parameter regime of the classical 
spin models. However, we remark that the entries 
of our entries correspond to the Boltzmann weights 
of the interactions, which are not only real but also 
positive. The latter condition is not satisfied in [34] . 

(ii) We observe that universality is already obtained 
for a suitable discrete set of unitary gates e*^^®^ 
[32, 33], which leads to a discrete set of Boltzmann 
weights in the corresponding six vertex model. 



B. Ising model 

Here we show that approximating the Ising partition 
function in the presence of an external magnetic field is 
a BQP-complete problem. More precisely, we find the 
following. 

Result 2 (BQP— completeness of the Ising model) 

Consider any planar circuit graph G. Let r denote the 
number of horizontal edges in G and let n be its vertical 
dimension. Consider a classical Ising model at inverse 
temperature P defined on G, where on each site a con- 
stant (complex) magnetic field ha is present satisfying 
^(3ha = and on each edge a constant (complex) cou- 



pling Je is present satisfying e 



Let Z denote the 



partition function of the model with open boundary con- 
ditions. Then we provide efficient quantum algorithms to 
estimate 



Z 



2^ 



+n 



(41) 



with polynomial accuracy. We also show that the problem 
of estimating (41) with this accuracy is BQP-complete. 

The proof will consist of several steps. Gates corre- 
sponding to this model are of the form (23); in particular, 
we consider the gates 



Wh := 
V := 



i 1 
1 i 



Wy := diag(z, 1, l,i) 



i7r/4 Q 

1 



(42) 



To show that (41) can be approximated with polyno- 
mial accuracy in poly-time with a quantum computer, 
we use the mapping of an edge model in the presence of 
a local magnetic field (with open boundary conditions) 
to a quantum circuit C described in Sec. IV B. Letting n 
be the vertical dimension of G as in the statement of the 
result, the associated quantum circuit C is an n-qubit cir- 
cuit composed of the gates Wh^ Wy and V. In particular, 
we have 



where |+) = (|0) + |1))/V^. Note that, for the values 
of the magnetic fields and the couplings adopted in the 
statement of the result, the matrices Wy and V are uni- 
tary. Moreover, the matrix Wh '= Wh/V^ is unitary as 
well. Letting C denote the unitary quantum circuit ob- 
tained by replacing every gate Wh by Wh, we simply have 
C = C/2^. It follows that 



Z 



2t 



+n 



(44) 



where the right hand side now represents a matrix ele- 
ment of a poly-size, unitary quantum circuit. From the 
Hadamard test it follows that estimating (41) with poly- 
nomial accuracy is achievable in poly-time on a quantum 
computer. 

Next we show that approximating (41) with polyno- 
mial accuracy is BQP-hard. This proceeds in the follow- 
ing steps: 

1. Denote T := V^WhVi, where Vi := diag(e^,l). 
We will need the following result: 

Lemma 1 (Universality of Ising— type gate set) 

Any poly-size n-qubit quantum circuit composed of two- 
qubit unitary gates can be approximated to accuracy e 
(with respect to the operator norm) by a circuit composed 
of poly{n, ^) single-qubit gates T and two-qubit gates 
T^'^WyT^'^ , where every such two-qubit gate is restricted 
to act on nearest-neighboring qubits only. Moreover the 
latter circuit can be found in time poly{n^ 

This Lemma is proved in Appendix A; here we con- 
tinue the argument to prove Result 2. 

2. Consider an arbitrary poly-size n-qubit circuit 
U composed of two-qubit unitary gates. It fol- 
lows from the discussion in Sec. Ill that the 
problem of approximating general matrix elements 
(+|^^/7|+)^^ with polynomial accuracy is BQP- 
hard. We now write U as 



U = [vifU'[V^] 



(g)n 



(45) 



for some suitable [/^ Note that is again a poly- 
size circuit, and can be found efficiently. We denote 
A:= [Vi]^^, i.e. U = AU'A. 

3. Due to Lemma 1, every poly-size circuit can be 
approximated with accuracy e by a circuit U'^ of 
size poly(n, 1/e) composed of T gates and nearest- 
neighbor T^'^WyT^'^ gates, and U^^ can be found 
efficiently. This implies that each matrix ele- 
ment (+|*^"^/7|+)*^^, where U is any poly-size cir- 
cuit as before, can be approximated with accuracy 
l/poly(n) by a matrix element of the form 



(46) 



Z = 2^(+|^^C|+) 



(43) 



where U^' is a poly-size circuit composed of T gates 
and nearest-neighbor T'^'^WyT^'^ gates. Hence, 
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also the problem of approximating the matrix ele- 
ments (46) with 1/poly accuracy is BQP-hard. Our 
goal is to show that any such matrix element coin- 
cides with a quantity Zj k^'^ considered in Result 2. 

4. By definition, the gates T and T®^W^T®^ contain 
square roots of the operation V. However, the op- 
eration AU" A only contains integral powers of the 
gate V . To see this, first note that any circuit U" 
of T and T®^W^T®^ gates only contains gates 
at the left and right 'boundary' of the circuit. In 
the circuit AU" A^ each V2 gate at the boundary is 
multiplied with another V2 gate due to the pres- 
ence of the left and right boundary operator A. As 
a result, the operation AU' A is a genuine circuit 
composed of W^, Wh and V gates. 

5. With any such circuit AU" A we now associate a 
graph G in the following way. With each single- 
qubit gate Wh we associate a single horizontal edge, 
and with each two-qubit gate we associate a 
vertical edge. The entire graph G associated with 
the quantum circuit is then obtained by simply 'glu- 
ing' together these edges in the natural way. An 
example is given in Fig. 4. Note that the graph is 
indeed a planar circuit graph, and that the vertical 
dimension of G is n. Moreover, consider an Ising 
model defined on G with parameters /3, Jg and /ig 
as in Result 2, with open boundary conditions, and 
let Z denote the partition function of the model. 
Then Z can be expressed as a matrix element of 
the form (14), for some circuit C composed of W^^ 

and V gates. It is now straightforward to ver- 
ify that, up to a normalization stemming from the 
presence of 1/a/2 in the gates W^^ the circuit C 
coincides with the circuit AU" A. More precisely, 
letting r denote the number of horizontal edges in 
G, one has C = 2^ AIJ" A. This shows that 



{-\-\®'' AU" A\ 



Z 



2f +n • 



(47) 



It follows that every matrix element of the form 
(46) coincides with an Ising partition function with 
parameters /3, Jg and ha as in Result 2, defined on 
the associated planar circuit graph G. This proves 
that the problem of estimating Z /2^^'^ with poly- 
nomial accuracy is BQP-hard. 



C. Potts model 

In this section we show that approximating partition 
function of the three-level Potts model with complex pa- 
rameters on a quasi 2D square lattice is BQP-complete. 
To be precise, we find the following. 

Result 3 (BQP— completeness of the Potts model) 

Consider the Potts model with three-level particles de- 
fined on a poly-size quasi 2D square lattice with fixed 



boundary conditions. Let ^ denote the inverse tempera- 
ture and Ju=v (Ju^v) the coupling strength when the two 
interacting particles spins u and v are (not) in the same 
state (see (25)^). Consider this model defined in the fol- 
lowing parameter regime 



(e 



{(l,0),(e-/«,l),(-i,l), 
(l,l),(0,l),(e,l), 

((V2e)-\l),(-l,l)} (48) 



where e is a polynomially small number, i.e. e = 
(9(l/poly(n)). Moreover, consider that certain 'blocks^ 
of these coupling strengths appear together (that is, there 
are certain local distributions of couplings). Let Z denote 
the partition function of this model. Then we provide effi- 
cient quantum algorithms to estimate Z with polynomial 
accuracy. We also show that estimating Z with this ac- 
curacy is BQP-complete. 

To prove this result, we will make use of the idea of 
encoded universal quantum computation as in Sec. VA. 
We will first present the encoding of physical qubits into 
logical ones, and then we will construct a universal gate 
set at the logical level with Potts-type gates. 

First of all we define an encoding. According to our 
mapping of Sec. IV B, a three-level quantum particle (a 
qutrit) is associated to each three-level classical particle. 
Let {|0), |1), |2)} denote the basis states of each qutrit. 
We use two of these qutrits (which we refer to as physical 
qutrits) to encode one logical qubit, whose states |0) and 
|1) are defined as 



ID 



|o>|i>, 

|1>|2>. 



(49) 



We will refer to the first and the second physical qutrits 
as the 'upper' and the 'lower qutrits', respectively. Note 
that, with this encoding, the input state |0)|0) . . . |0) re- 
quires no preparation, since it can be fixed as an initial 
boundary condition, viz. = |0) |1) . . . |0) |1). 

Now we proceed with the construction of the encoded 
universal gate set. We first observe that, due to the Potts 
interaction (25), the gates (16) and (17) contain only 
two different coefficients; for further reference, we denote 
them 

:= (e^J-^^e^^-^-). (50) 
We construct the encoded universal gate set as follows: 



• The single-qubit identity 

1 

ii= Yl , 

i,j=0 



(51) 



is achieved by applying a two-qutrit gate (between 
the two physical qutrits that compose the logical 
qubit) with (/i, u) = (1, 1) (see Fig. 6). 
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• The phase gate 

P^/8 = |0){0|+e-/«|l){l| (52) 

is achieved by applying a two-qutrit gate between 
an auxiliary qutrit in the state |2) and the lower 
qutrit of the logical qubit with (/i, z^) = (e*^/^, 1) 
(see Fig. 6). 

• The Hadamard gate 

^ = 4 E (-ir^HH (53) 

^ u,v=0 

is more involved and we refer the reader to Ap- 
pendix B for the details. 

• The two-qubit identity gate 

1 

I2 = Yl (54) 

u,v=0 

is trivially obtained by applying a two qutr it-gate 
between the two physical qutrits of the first logical 
qubit with (/i, z^) = (1,1), and doing the same for 
the other logical qubit (see Fig. 7). 

• Finally, the phase gate between logical qubits 

1 

CZ= Y \uv){uv\ (55) 

u,v—0 

is achieved by applying a two qutrit-gate between 
the lower qutrit of the first logical qubit and the up- 
per qutrit of the second logical qudit with (/i, u) = 
(—1, 1). We also need to apply the same gate be- 
tween the lower qudit of the second logical qubit 
and an auxiliary particle in the state |2) (see Fig. 7). 




(a) (b) 



FIG. 6: Auxiliary qubits (i.e. fixed qubits) are depicted in 
gray, physical qubits in black, and logical qubits are high- 
lighted in yellow, (a) Logical single-qubit identity Ii. (b) 
Logical phase gate P-n/s- Each gate is determined by the pair 
of numbers (/x, v) (Eq. (50)) which are indicated next to it, 
and its color is just a guide to the eye to identify equal gates. 
Note that time runs in both figures from right to left. 

In the construction presented above we need a constant 
supply of auxiliary qutrits, for which we envisage the 
following setup. We arrange the particles so that each 



(1,1) 




^ ^(-1,1 ) 

time time 
(a) (b) 

FIG. 7: (a) Logical two-qubit identity I2. (b) Logical con- 
trolled phase gate CZ. See the caption of Fig. 6 for the 
meaning of the pair of numbers and the color associated to 
each gate. 



physical qubit is in contact with four auxiliary qubits, 
as indicated in Fig. 8. Moreover, each physical qubit is 
propagated in time via a single-qubit gate at the physical 
level (that is, a gate of the form (16)) which is achieved 
with {ji^f) = (1,0)- In this way, the physical qubits are 
distributed on a 2D square lattice, which is surrounded 
(in front and on the back) by auxiliary qubits — we refer 
to this as a quasi 2D square lattice. 

This proves that we can perform encoded univer- 
sal quantum computation using qutrits interacting with 
Potts-type nearest-neighbor gates. From the discussion 
of Sec. Ill it follows that estimating (0|^^C|0)^^ = Z 
with polynomial accuracy is BQP-complete, where Z is 
the partition function specified in the statement of Result 
3. This thus concludes the proof of Result 3. 

We point out that our BQP-completeness result for the 
Potts model imposes constraints in the parameter regime 
as noted in Sec. IV B. We now give some further remarks 
regarding this proof. 

(i) Note that to realize a gate in more than one step 
(like the Hadamard gate, see Appendix B) implies 
that one must fix the value of all those coupling 
strengths together in order to realize that gate in the 
quantum circuit. That is, the gate is only realized 
(and thus the complexity result is only achieved) if 
the classical model contains that 'block' of coupling 
strengths. These are thus constraints in the distri- 
bution of couplings imposed by the realization of 
gates. 

(ii) We observe that the single-qubit identity at the 
physical level /i requires to set = 0, which for- 
mally corresponds to setting Ju^v = —00. This 
unphysical regime can be avoided by letting e^'^"^^ 
be a polynomially small number, that is e^'^"^^ = 
(9(l/poly(n)) By the same argument as for the 
Hadamard gate (see Appendix B), the resulting 
state will be polynomially close to the ideal one, 
and thus the error within the accuracy of the com- 
putation. 
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phys. 




|2> 



|0> 



(1,1) 



(1,0) 



(a) 



phys 

aux. I 



aux. 





FIG. 8: Setup for a Potts model with the physical qutrits 
surrounded by auxiliary qutrits (i.e. qutrits whose value is 
fixed). Single-qutrit gates are applied in the time direction, 
and two-qutrit gates are applied in a fixed time slice. Solid 
and dashed lines stand for the two-qubit and single-qubit 
identities (at the physical level), respectively, as specified in 
the legend, (a) A time slice of the lattice at time t, where one 
can see the structure with two auxiliary qutrits connected to 
each physical qutrit. Each logical qutrit is composed of two 
physical qutrits. (b) Two time slices of the lattice, at times t 
and t + 1. The figure shows an example of a circuit, where a 
phase gate Ptt/s is applied to and the gate CZ is applied 

to 



(iii) We see with our construction that the Potts model 
in a ID array is not likely to be BQP-complete, 
whereas this model on a 2D setup (and with com- 
plex parameters) is BQP-complete. This is in agree- 
ment with the fact that the Potts model in ID array 
or a tree-like structure is efficiently simulable clas- 
sically [35]. 

(iv) Finally, this result can be generalized to Potts 
model with any q. In this case, the encoded states 
are still given by Eq. (49), and all gates are per- 
formed with the same procedure except for the 
Hadamard. For this gate, one adds filters for the 
\3) . . . \q — l) components on the upper and the lower 



qudit right after the filters for the |2) and the |0) 
component of Fig. 10. That is, the number of filters 
scales linearly with q (since one requires q — 2 filters 
for each qudit). Moreover, each physical qudit has 
to be connected to 2\q/2] auxiliary particles, each 
fixed in a different state, namely |0),...,|g — 1). 
This amounts to a similar construction to that of 
Fig. 8 but where each physical qudit is connected to 
\q/2] auxiliary qudits to the left and to the right. 



D. Z2 Lattice gauge theory 

Now we turn to a different class of models, namely 
LGTs. Using the tools presented in Sec. IV C we will 
prove the following result. 

Result 4 (BQP-completeness of the 3D Z2 LGT) 

Consider a Z2 lattice gauge theory defined on a 3D rect- 
angular lattice of size (4x, 12y^ 7z), where x, y and z are 
natural numbers, and with fixed boundary conditions. Let 
(3 denote the inverse temperature and Jf the coupling 
strength on its faces. Consider this model in the following 
complex parameter regime, 



(56) 



where is a continuous parameter ^ G [0,27r); and ( is 
a polynomially small number, ( = 0(1/ polyin)). More- 
over, consider that certain 'blocks^ of couplings appear 
together, and that there is a certain gauge fixing. Let 
Z denote the partition function of this model. Then we 
provide efficient quantum algorithms to estimate 



Z 



2^ 



(57) 



with polynomial accuracy. We also show that the problem 
of approximating (57) is BQP-complete. 

The proof of this result will proceed similarly as in the 
previous sections: we will construct an encoded universal 
gate set. We first define the following encoding of four 
physical qubits into one logical qubit: 



|0) = |0)|0)|0)|0) 

ll> = ll)ll>ll>ll)- 



(58) 



Note that preparing the input state |0) is trivially 
achieved by fixing the physical qubits (which correspond 
to physical spins) in |0) |0) |0) |0). In the following we show 
how to construct a general single-qubit unitary gate (at 
the logical level), and a specific two-qubit gate (also at 
the logical level), namely (the non-local part of) a con- 
trolled phase gate. 

• First we express an arbitrary single-qubit rotation 
in its Euler decomposition. 



t7(7,/3,a) = R,{^)HR,{l3)HR,{a). 



(59) 
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Thus, our goal is to generate arbitrary z-rotations 
and the Hadamard gate. The gate Rz{£,) can be 
implemented by letting the logical qubit interact 
with a neighboring, auxiliary face which has all re- 
maining qubits fixed to |0), and setting e^'^ = e*^ 
in this face (see Fig. 9). 

The logical Hadamard gate H is implemented as 
a teleport at ion-based gate. We refer the reader to 
Appendix C for the details. We shall in fact apply a 
non-normalized version of this, namely H := ^/2H 
(see below for a discussion of this normalization 
factor). 

• To generate the logical controlled phase gate CZ 
we decompose it into single-qubit Rz rotations and 
a two-qubit gate: 

CZ = (- 1 (- 1) diag(l, (60) 

The two-qubit gate diag(l, i, i, 1) is implemented 
by letting the last two logical qubits interact via an 
auxiliary plaquette with e^^ = i (see Fig. 9). 




FIG. 9: The physical qubits composing the logical qubit are 
marked with thick, black lines, and wavy lines indicate qubits 
fixed by the gauge symmetry. Each gate is determined by 
e^'^. (a) The logical rotation around the ^-axis, Rz{0 
obtained by letting the logical qubit interact with an auxiliary 
plaquette whose spins are all fixed to |0), and setting e^'^ — e*^ 
in this auxiliary plaquette. (b) The non-local part of CZ, 
diag(l, i, i, 1), is implemented by letting the two logical qubits 
interact via an auxiliary plaquette with e' 



(3J 



Finally, we apply single-qubit identities at the physical 
level Ii := |0)(0| + |1)(1| on all physical qubits which be- 
long to a logical qubit, unless otherwise specified. This 
gate is achieved by rescaling the interaction from e^'^f 
to e^'^'^f~^\ and then letting e^*^-^^"^^ tend to polyno- 
mially fast (see Appendix C). We also set J = on the 
faces whose interaction is not mapped to a quantum gate, 
which sets their corresponding Boltzmann weights to 1, 
thus not 'contributing' to the partition function. This is 
in accordance with the fact that these particles do not 
take part in the computation described by the quantum 
circuit. 



Our construction of a circuit at the logical level dif- 
fers from the general mapping presented in Sec. IV C in 
the following sense. The quantum circuit processes a ID 
array of logical qubits which are distributed along the x 
direction. These are transformed by four-qubit physical 
gates in the x direction (leading to logical single- and 
two-qubit gates), and by single-qubit physical gates in 
the z direction (which corresponds to time). The addi- 
tional dimension {y direction) is required for the realiza- 
tion of the teleportation-based Hadamard gates, which 
transport the ID array that stores the quantum informa- 
tion in the y direction. 

Thus, we have constructed an encoded universal quan- 
tum circuit C with Z2-LGT-type gates, where C contains 
the non-normalized version of the Hadamard gates H. 
Let C denote the unitary circuit containing the normal- 
ized Hadamard gates H = H / \/2. Then we have that 



Z 



2xyz/2 



(0|^^C|0) 



(61) 



where xyz is the number of times that the gate H is 
applied. From (61) and the discussion of Sec. Ill it fol- 
lows that the estimating the partition function (57) with 
polynomial accuracy is BQP-complete. This concludes 
the proof of Result 4. 

We emphasize that, as in the previous results, we re- 
quire complex parameters due to the unitary gates. Also 
the remark (i) of Sec. VC applies here: performing a 
logical gate in several steps (like the Hadamard gate, see 
Appendix C) implies that a certain distribution of cou- 
plings must appear together in the classical model in or- 
der to prove our complexity result. Next we give some 
more comments concerning this result. 

(i) Note that the single- and two-qubit identity 
gates Ji and I2 can be generated using the 
universal gate presented above. For example, 
Rz{0)HRz{0)HRz{0) = Ii and = I2. 

(ii) Result 4 is in contrast with the 2D Z2 LOT, whose 
complexity is 'trivial', i.e. it is in P when the tem- 
poral gauge is fixed. On the other hand, the 3D Z2 
LGT can be mapped via a duality transformation 
to the 3D Ising model [20]. Thus, they share the 
same complexity in the parameter regime given by 
this transformation. We also observe that it was 
recently shown that computing the partition func- 
tion of the 4D Z2 LGT in a real parameter regime 
is #P-complete [36, 37]. Concerning the quantum 
computational complexity of the 4D Z2 LGT (and, 
more generally, of any d Z2 LGT with d > 3), one 
can make the following trivial observation: any of 
these models can be "reduced" (i.e. become effec- 
tively equivalent) to the 3D Z2 LGT if the couplings 
in every face except those in a 3D volume are set 
to 0. Thus, our results imply that their complexity 
in this (trivial) parameter regime is BQP-complete. 
However, to the best of our knowledge, no quantum 
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computational results are known outside this trivial 
parameter regime. 



VI. FURTHER RESULTS: ONE CLEAN QUBIT 
MODEL 

Here we point out a connection between the results 
obtained on the Ising model in Sec. VB and a scheme 
for quantum computation called the 'one clean qubit 
model' [25]. In the latter, one considers a quantum 
computation where all qubits but the first one are ini- 
tialized in the totally mixed state (and the first qubit 
is, say, in the state |0)). To this initial state, an arbi- 
trary poly-size quantum circuit may be applied, followed 
by a single-qubit measurement in the final stage of the 
computation. This one-clean-qubit model comprises a 
scheme that is believed to be weaker than the full power 
of quantum computers but stronger than classical com- 
putation (although these are unproved assertions). The 
corresponding complexity class of decision problems that 
can be solved efficiently with the one-clean-qubit scheme 
is called DQCl. 

A standard problem that can be solved using the one 
clean qubit model is the problem of estimating normal- 
ized traces of unitary quantum circuits. Let U denote 
a poly-size n-qubit quantum circuit composed of, say, 
two-qubit gates. Then there exists an efficient quantum 
algorithm within the one clean qubit scheme which re- 
turns a number c that provides (with exponentially small 
probability of failure) an e-approximation of the normal- 
ized trace 2~^Tt{U) in poly(n) time, for every e that 
scales at most inverse polynomially with n. The tech- 
nique is a simple variant of the Hadamard test and we 
refer to the literature [25]. Moreover, the problem of es- 
timating normalized traces of unitary quantum circuits 
is known to be DQCl-hard. In other words, this problem 
captures all problems that can be solved efficiently within 
the one clean qubit paradigm. It thus plays a similar role 
as the unitary matrix element problem for BQP. 

Our results obtained in Sec. VB immediately lead to 
a complete problem for DQCl involving Ising partition 
functions defined on graphs with periodic boundary con- 
ditions. We limit ourselves to a sketch of the argument. 
It follows from the discussion in Sec. VB that, for every 
poly-size quantum circuit U there exists an Ising model 
defined on a planar circuit graph G (which can be found 
efficiently), with couplings and temperature as in Re- 
sult 2 and with open boundary conditions, such that Z/k, 
provides a 1/poly approximation of the matrix element 
(+|^^/7|+)^^. Now, instead of |+)^^ we consider the 
matrix element {s\U\s) where \s) = |5i . . . Sn) represents 
an arbitrary computational basis state. This matrix el- 
ement now coincides with j where Z^ denotes the 
partition function of the same model (i.e. same graph, 
couplings and temperature) , however considering bound- 
ary conditions where the left and right boundary spins 
are fixed in the same configuration s. Such a situation 



is equivalent to considering a graph G' where each left 
boundary spin at the /cth 'row' in the graph is identified 
with its corresponding right boundary spin, and the re- 
sulting spin is fixed in the state s/^, for all k. One thus 
arrives at an Ising model on a new graph G' which is 
obtained from G by enforcing periodic boundary condi- 
tions^ and where one vertical slice of spins is fixed in the 
configuration s. Finally, consider the normalized trace 
of the matrix element [/, given by 2~^ Due 
to the above discussion, this normalized trace may be 
approximated with 1/poly accuracy by 

s 

But as Z^ is the partition function on G' with one ver- 
tical slice of spins fixed in the configuration x, the sum 
Z^ = Z' simply represents the partition function on 
G' where these spins are now fully summed out. The 
quantity Z' is hence the full-fledged partition function 
on G' of the Ising model with couplings and tempera- 
ture as before, and without any fixed spins or boundary 
conditions. 

We thus arrive at the following result: consider any 
planar circuit graph G. Let r denote the number of 
horizontal edges in G and let n be its vertical dimen- 
sion. Let G' be the graph obtained by enforcing periodic 
boundary conditions on G (as above). Consider a classi- 
cal Ising model at inverse temperature /3 defined on G, 
where on each site a constant (complex) magnetic field ha 
is present satisfying e^^^ = e^, and on each edge a con- 
stant (complex) coupling Jg is present satisfying e^^^ = z. 
Let Z' denote the partition function of the model. Then 
the problem of approximating 



with polynomial accuracy, is DQCl-hard. 

VII. CONCLUSIONS 

In this work we have shown that estimating the par- 
tition function of the six vertex model on a 2D square 
lattice, the 2D Ising model with magnetic fields on a pla- 
nar graph, the Potts model on a quasi 2D square lattice, 
and the 3D Z2 LOT is BQP-complete. Ah these mod- 
els must be defined in a (partially) complex parameter 
regime in order to prove the result. In all but the Ising 
model, we further require that certain blocks of coupling 
strengths appear together. Roughly speaking, that these 
problems are BQP-complete means that they are as hard 
as simulating arbitrary quantum computation. Because 
our proof is constructive, we have also provided the ef- 
ficient quantum algorithms that estimate the partition 
functions of the above models. The proofs are based on a 
mapping from partition functions to quantum circuits in- 
troduced in [23] and extended here to the standard Ising 
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model, Potts model and to Z2 LGTs. In this sense, our 
work puts these different kinds of classical spin model on 
an equal footing as far as their quantum computational 
complexity is concerned. 

It would be interesting to obtain quantum algorithms 
for computing such partition functions in a real param- 
eter regime. At present we do not know whether our 
approach can be extended to prove such results. As we 
pointed in Sec. VA remark (i), the Boltzmann weights 
give rise not only to real but also to positive entries in the 
quantum gates. The latter seems to be a rather severe 
restriction in the pursuit of this goal. 
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Appendix A: Universality of the Ising— type gate set 



The proof of Lemma 1 will use the following lemma. 



Lemma 2 Let e > and define K := WhV . Up to a 
global phase, any two-qubit unitary operation can be ap- 
proximated with accuracy e (with respect to the operator 
norm) by a product of gates K and Wy i.e. these gates 
form a strongly universal gate set. 

Proof: Let P := diag(l, i) denote the phase gate, H the 
usual Hadamard gate, H = j=o ^ 
denote the standard Pauli gate. We will also denote 
the controlled phase gate CZ := diag(l, 1, 1, — 1). Then 
the following identities can easily be verified: 



ZKZKZ{K'^)^Z (X H 
KZK^ZK^ZK 



(Al) 



Now note that (x Z ® Z. Thus Eqs. (Al) imply that 



(A2) 



So far we have showed that the operations H ^ I and 
P I can be written as a product of gates K, and 
Wy] the same is easily seen to hold for / (g) and 1 P. 
Consequently: 

(i) As CZ oc [P0P]Wv, the CZ gate can also be written 
as a product of and Wy gates; 



(ii) As Wl (X HP^H, the gate V 



(regarded as 



a two-qubit operation acting nontrivially on either 
the first or second qubit) can be written as a product 
of and Wy gates as well; 



(iii) As the gate K does not have finite order, there does 
not exist any integer m such that = . How- 
ever, it is well known that for every ^ > there 
exists an integer m = poly(l/(5) such that the dis- 
tance (in operator norm) between and is at 
most S. 

This shows that the gates i7, U as well as the CZ gate 
can be approximated with accuracy e by a product of 
poly(l/e) gates K and Wy. Moreover, it is known that 
the CZ, H and V gates generate a dense subgroup of 
/7(4) (up to global phases). This hence proves Lemma 2. 
□ 

The proof of Lemma 1 is now obtained as follows. We 
simultaneously conjugate the gates K and Wy with the 
operation acting on each qubit, yielding 



V2 U2 



Wy[Vi 0U^]t 



T 

Wy. 



(A3) 



As K and Wy generate a dense subgroup of /7(4) up to 
global phases due to Lemma 2, the same holds for T and 
Wy since the latter gates are obtained by conjugating 
the former with the same fixed local unitary operation. 
Finally, it follows that also the gates T and T^'^WyT^'^ 
generate a dense subgroup of [/(4) up to global phases: 
indeed, can be approximated with arbitrary accuracy 
by powers of T, such that Wy = T^^ [T^^WyT^^jT^^ can 
be approximated with arbitrary accuracy by products of 
T and T^'^WyT^'^. 

We have thus showed that every two-qubit operation 
can be approximated with arbitrary accuracy by prod- 
ucts of T and T^'^WyT^'^. The Solovay-Kitaev Theo- 
rem then guarantees fast convergence, i.e. every poly-size 
n-qubit circuit composed of two-qubit gates can be ap- 
proximated with accuracy e by a circuit of poly(l/e,n) 
gates T and T'^'^WyT^'^; moreover, the latter circuit can 
be found efficiently, i.e. in poly(l/e, n) time. Hence this 
proves the claim. 



Appendix B: Hadamard rotation with Potts— type 
gates 

In the following we show how to obtain a Hadamard 
gate with Potts-type gates. Consider a general state (de- 
fined at the logical level), = a|0) + where a and 
/3 are arbitrary, normalized coefficients, |ap + = 1. 
Our goal is to apply a Hadamard gate to this state, that 
is, to transform it to := + — ), where 



1 

71 



|0>±|1)). 



(Bl) 



by analogy with the conventional definitions |±) := (|0)± 
\1))/V2. 

To apply this gate, we first observe that a Hadamard 
rotation can be decomposed as 



H = R^{3TT/2)Ra,{TT/4)R^{3Tr/2) . 



(B2) 
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Then these three gates can be achieved with the fohowing 
sequence of Potts-type gates (see Fig. 10): 

1. A two qutrit-gate between an auxihary qutrit in 
the state |2) and the lower qutrit of the logical 
qubit with (/i, i^) = (— i,l), which results in the 
state a|0)|l) -ip\l)\2). 

2. A single-qutrit gate on the upper qutrit with 
{(JL^f) = (— 1), and then a single-qutrit gate on 
the second qutrit with (/i, u) = (1,1), which results 
in the state a{-i\0) + |1) + |2))(|0) + |1) + |2)) - 
z/3(|0)-z|l) + |2))(|0) + |l) + |2)). 

3. A two-qutrit gate between an auxiliary qutrit in 
|2) and the upper qutrit with (/i, z^) = (0,1), and 
the same gate between an auxiliary qutrit in |1) and 
the lower qubit, which results in the state a{—i\0)-\- 
|l))(|l) + |2))-i/3(|0)-i|l))(|l) + |2)). 

4. A two qutrit gate between the upper and lower 
qutrit with (fi^v) = (0,1), which yields the state 
a(-i|0)|l) - i|0)|2) + |1)|2)) - i/?(|0)|l) + |0)|2) - 
i|l)|2)). 

5. A two qutrit gate between an auxiliary qutrit in 
|0) and the upper qutrit with (/i, z^) = (e, 1), and 
the same gate between an auxiliary qutrit in |2) 
and the lower qutrit. Then a two-qutrit gate be- 
tween an auxiliary qutrit in |1) and the upper qutrit 
with {jJ^jv) = (1/(6^/2), 1), and the same gate be- 
tween an auxiliary qutrit in |1) and the lower qutrit. 
Up to terms of order O(e^), the resulting state is 

[aH|0)|l) + |1)|2)) - i/3{\0m - i\m)]/^- See 
below for a discussion of the accuracy of the gate 
due to neglecting these higher order terms. 

6. A two-qutrit gate between an auxiliary qutrit in |0) 
and the upper qutrit with (^, 1), and another two- 
qutrit gate between an auxiliary qutrit in |1) and 
the upper qutrit. The resulting state is [a(|0)|l) + 
i|l)|2)) + /?(|0)|l)-i|l)|2))]/V2. 

7. Finally, a two-qutrit gate between an auxiliary 
qutrit in 1 2) and the lower qutrit with (— 1), after 
which the state of the system is [a(|0)|l) + |1)|2)) + 
/3(|0)|1) - |l)|2))]/v/2 = a{\0) + \1))/V2 + /3(|0) - 
|1))/a/2 = H\%l:)^ as we wanted to show. 

Note that, due to step 5, the desired gate is achieved up 
to an error of C^(e^), where we use e.g. the Jamiotkowski 
Fidelity F [38] as a measure for the gate fidelity. Mak- 
ing use of the distance measure between the imper- 
fect operation and the ideal Hadamard gate, D{E,U) = 
Vl — together with the chaining inequality [38, 39], 
D{£io82Mi 0^2) < D{8iMi)^D{82M2), we conclude 
that the total distance after M applications of an imper- 
fect operation of this kind in a quantum circuit is at most 
M times the distance of a single imperfect operation. 
That is, \/l — Ftot < — F, where Ftot denotes the 



Jamiolkowski fidelity of the overall circuit containing M 
imperfect operations, each with Jamiolkowski fidelity F. 
Hence, for any circuit consisting of a polynomial number 
of gates, M = 0{N)^ a polynomial accuracy in e, that is 
e = 0(l/poly(A^)), suffices to achieve a polynomial ac- 
curacy of the total circuit. Notice that this is sufficient 
to obtain a final approximation of the partition function 
with polynomial accuracy [42]. 
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FIG. 10: The Hadamard rotation at the logical level, H, 
is achieved by first performing a z-rotation Rz, then an x- 
rotation Rx (which is itself composed of various gates) and 
finally another z-rotation Rz . Most of the gates involve aux- 
iliary qubits, whose input state is indicated on the right (since 
time runs from right to left). See the caption of Fig. 6 for the 
meaning of the pair of numbers and the color associated to 
each gate. See the text for the transformation of the state |'0) 
at each stage of the circuit. 



Appendix C: Hadamard rotation with Z2 Lattice 
Gauge Theory— type gates 

Here we show how to general single-qubit rotation of 
(59) with Z2 LGT-type gates. This operation includes 
Rz rotations, which we have shown how to perform in 
Sec. VD, and the non-normalized version of the logical 
Hadamard gate H. We will implement the latter as a 
teleportation-based gate. Intuitively, the idea is the fol- 
lowing. To apply H to the logical qubit l^^), we prepare 
two logical qubits next to it in the state |0+) + |1 — ) 
(see (Bl)). Then we apply a projection onto a Bell state 
1 00) + 111) between the qubit we want to teleport and the 
first of these two logical qubits. The third logical qubit 
is then left in the desired state, namely H\il)). 

More precisely, consider a general initial logical state 
= a|0) + where a and h are unknown, normal- 
ized coefficients |ap + \h\^ = 1. Then we transform this 
state into Rz{7r/2)HRz{a)\ijj) by applying the following 
sequence of steps (see Fig. 11): 



1. Single-qubit gates with 



1 to all physical 



qubits of the two logical qubits on the right. Each 
of these gates transforms the state of the physical 
qubits from |0) to |0) + |1) =: V2|+). Thus, this 
results in the state |'0)16| + + + +)^^. 
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2. The four-qubit gate with e^^ = e^^^+^Z^^ on the 
top physical qubit of the first logical qubit. Also the 
four-qubit gate with e^^ = inside the plaquette 
where the last two logical qubits are defined, and 
the four-qubit gates with e^^ = on the physical 
qubits on top and on the right of the last two logical 
qubits (as indicated in the second step of Fig. 11). 
This results in the state (a|0)+e^^"+^)6|l))[|0000) + 
|0011) + |1100) + |1111)]^2^ 

3. The four-qubit gates with e^^ = between the 

physical qubits on the right and on the bottom 

of the last two logical qubits (as indicated in the 

third step of Fig. 11). The resulting state is 
(a|0)+e^("+")6|l))[|0) + |l)]^2 

4. The four-qubit gate with e^^ = i applied on the 
top physical qubit of the second logical qubit. A 
four-qubit gate with e^^ = between the physical 
qubit on the right of the second logical qubit and 
the physical on the right of it, which we refer to 
as the "middle qubit" (see fourth step of Fig. 11), 
which propagates the state of the first qubit to the 
middle one. Then the four-qubit gate with e^'^ = 
i between the middle qubit and the left physical 
qubit of the third logical qubit. The resulting state 
is (a|0)+e^(-+-)6|l))[|0)(|0)+z|l))-|l)(|0)-z|l)]. 

5. Finally, two concatenated four-qubit gates with 
g/3j _ Q ]3etween the first and the second logical 
qubit, which corresponds to the projection onto 
the Beh state (0|(0| + (1|(1|. We are interested 
in the state of the third logical qubit, since this 
contains the result of the teleportation. This state 
is a(|0) + + e^"6(|0) - z|l)), which equals 

R,{^/2)HR,{am. (CI) 

Thus, this sequence of gates applies essentially the 
first two logical gates of a general single qubit ro- 
tation of (59). The rotation Rz{7t/2) appearing in 
(CI) can be compensated in the next rotation (marked 
as time step 8 in Fig. 12) by applying a four- 
qubit gate with e*^^+^/^^ (instead of e*^^+^\ which 



would be the analogous case to the rotation of step 
3, where 6**^^+^^ was applied). This results in the 
state Rz{P)HRz{(^)\'^)- Then a teleportation-based 
Hadamard gate H is applied with the same procedure as 
the one described above, and this gives rise to the state 
Rzii^ /2)HRz{l3)HRz{a). The final rotation must cor- 
rect for this phase, and the four-qubit gate at step 11 of 
Fig. 12 must have e**^^~^/^^ (instead of e^^). The over- 
all sequence implements the desired general single qubit 
gate, R,{j)HR,{l3)HRMM)- 

Moreover, we apply single-qubit identities at the phys- 
ical level Ii := |0)(0| + |1)(1| on all physical qubits 
which belong to a logical qubit unless otherwise speci- 
fied, as mentioned in Sec. VD (yellow faces in the tem- 
poral direction in Fig. 11). This gate is achieved by 
first rescaling the face interaction energy (Eq. (5)) to 
— {Jf — l)S{si + Sj + 5/e + Si). According to (28) this 
gives rise to single qubit gate whose diagonal elements 
are 1 and the off-diagonal elements are e^^'^f~^\ Then 
the gate is achieved by setting e^^'^f~^^ = where ^ is 
a polynomially small number. By the same argument as 
for the Potts single-qubit identity (see Sec. VC remark 
(ii)), this will result in polynomial error at the end of the 
computation, which is within the accuracy of the result. 

Note that the unnormalized gate H is applied once in 
every volume of 4 units in the x direction, 12 units in the 
y direction, and 7 units in the z direction. As mentioned 
in Sec. VD, the number of these blocks is what gives rise 
to the factor appearing in (57). 

Since the general single qubit unitary gate involves a 
number of spins fixed by the gauge, it is important to 
verify that no closed loops of fixed spins are formed [28] . 
This is because the gauge fixing only yields an equivalent 
theory if no loops are formed, since, e.g., order param- 
eters of the theory are defined as products of variables 
on such closed loops. To verify that no such loops are 
formed we project our 3D lattice into a 2D lattice along 
the time direction (see Fig. 12), where one can easily rec- 
ognize that there is no closed path of edges fixed by the 
gauge. 
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FIG. 11: Part of a circuit with Z2-LGT-type gates that shows the processing of a logical qubit First, a logical rotation 
Rz{oi) is applied to it. Then, a logical teleportation-based Hadamard gate (up to a phase) is applied, which results in the state 
Rz{n/2)HRz{a)\'il^) (see text). Note that the qubit is processed in the z direction (corresponding to time) as well as the y 
direction. Colored faces indicate faces where the coupling strength has a specific value, as indicated in the legend. Wavy lines 
indicate edges whose spins have been fixed by the gauge, and their color has no meaning nor a correlation with the color of the 
faces. They are a guide to the eye to verify that no loops are formed, as shown in Fig. 12. 
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A planar graph is a graph which can be drawn in the 

plane without crossings of the edges. 

This is opposite to the way quantum circuits are usually 

drawn, where the first processing takes place from left to 

right. 

One may also use results from fault-tolerant quantum 
computation, which state that also noisy gates with a 
small but constant error suffice to perform quantum com- 
putation. This implies that e can be a small constant not 
depending on the size of the circuit, however additional 
overhead for encoding and fault- tolerant gate implemen- 
tation would be required in this case. 
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Two-qubit gate 



General single-qubit gate Rz{j)HRz{P)HRz{cx) 
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FIG. 12: A projection of the circuit of the 3D Z2 LGT into its spatial dimension. Part of the circuit in its normal form (that 
is, in the structure two-qubit gate, single-qubit gate, two-qubit gate) is illustrated. Logical qubits are indicated with thick, 
black lines. All edges in the time direction (going out of the paper) which are boundary to a logical qubit are fixed by the gauge 
(pink dots). Edges in the spatial direction are either not fixed by the gauge (black, thin lines), or are fixed by the gauge at 
different time steps 1 ... 13 indicated with a wavy, colored edges. Hence, a loop of edges fixed by the gauge would correspond 
in this figure to a closed loop of wavy, colored edges. It can be verified by inspection that no such loops are present. 



